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We show that, in principle, a slowly evolving gravitationally collapsing perfect fluid cloud can asymptotically 
settle to a static spherically symmetric equilibrium configuration with a naked singularity at the center. We 
consider one such asymptotic final configuration with a finite outer radius, and construct a toy model in which 
it is matched to a Schwarzschild exterior geometry. We examine the properties of circular orbits in this model. 
We then investigate observational signatures of a thermal accretion disk in this spacetime, comparing them with 
the signatures expected for a disk around a black hole of the same mass. Several notable differences emerge. 
A disk around the naked singularity is much more luminous than one around an equivalent black hole. Also, 
the disk around the naked singularity has a spectrum with a high frequency power law segment that carries a 
major fraction of the total luminosity. Thus, at least some naked singularities can, in principle, be distinguished 
observationally from black holes of the same mass. We discuss possible implications of these results. 
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I. INTRODUCTION 

The no hair theorems in black hole physics state that a black hole is completely characterized by three basic parameters, 
namely, its mass, angular momentum, and charge. Because of this extreme simplicity, the observational properties of a black 
hole are determined uniquely by these three intrinsic parameters plus a few details of the surrounding environment. 

There is compelling observational evidence that many compact objects exist in the universe, and there are indications that 
some of these might have horizons |1]. However, there is as yet no direct proof that any of these objects is necessarily a black 
hole. The dynamical equations describing collapse in general relativity do not imply that the final endstate of gravitational 
collapse of a massive matter cloud has to be a black hole (see |2]); other possibilities are also allowed. In the case of a continual 
collapse, general relativity predicts that a spacetime singularity must form as the endstate of collapse. However, recent collapse 
studies show that, depending on the initial conditions from which the collapse evolves, trapped surfaces form either early or late 
during the collapse. Correspondingly, the final singularity may either be covered, giving a black hole, or may be visible as a 
naked singularity. If astrophysical objects of the latter kind hypothetic ally form in nature, it is important to be able to distinguish 
them from their black hole counterparts through observational signatures. 

General relativity has never been tested in the very strong field regime, and very little is known about how matter behaves 
towards the end of the gravitational collapse of a massive star, when extremely large densities are reached and quantum effects 
possibly become relevant. Over the years, it has often been suggested that some exotic stable state of matter might occur below 
the neutron degeneracy threshold, allowing for the existence of quark stars, or boson stars, or even more exotic astrophysical 
objects (see for example yl). 

At present we cannot rule out the possibility that such compact objects do exist. On the other hand, the densities and sizes of 
compact objects in the universe vary enormously, depending on their mass. While a stellar mass quark star could have a density 
greater than that of a neutron star, a supermassive compact object at the center of a galaxy might have a density comparable to 
that of ordinary terrestrial matter. This means that it is difficult to come up with a single paradigm for all compact objects by 
simply modifying the equation of state of matter. Also, the processes that lead to the formation of stellar mass compact objects 
are different from those leading to supermassive objects. For one thing, time scales are very different. Stellar mass compact 
objects form in a matter of seconds, when the core of a massive dying star implodes under its own gravity. We know little about 
the physical processes that lead to the formation of supermassive compact objects, but whatever it is, it doubtless operates far 
more slowly than in the stellar-mass case. 

In our view it is important to study viable theoretical models that, under reasonable physical conditions, lead to the formation 
of different kinds of compact objects, and to investigate the properties of these different end states. In recent times, much 
attention has been devoted to the observational properties of spacetimes that describe very compact objects or singularities 
where no horizon is present ^. Many valid solutions of Einstein field equations exist which describe spacetimes that are not 
black holes. These are either vacuum solutions with naked singularities or collapse solutions in the presence of matter. Naked 
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singularity models include the Reissner-Nordstrom, Kerr and Kerr-Newman geometries with a range of parameter values that 
differentiate the black hole and naked singularity regimes. Collapse models include dust collapse, models with perfect fluids and 
those with other equations of state |2]. 

Note further that the classical event horizon structure of the Kerr metric can be altered in many ways. One way is by over- 
spinning the Kerr black hole in order to obtain a naked singularity [5J. Another is by introducing deformations or scalar fields 
to alter the spacetime and thus expose a naked singularity ^. Although the physical viability of some of these examples is not 
clear, the fact remains that classical general relativity allows for the formation of naked singularities in a variety of situations. 
If singular solutions represent a breakdown of the theory in the regime of strong gravity, then a study of some of these models 
might provide clues to new physics. For example, string theory or quantum gravity corrections can remove the Kerr singularity, 
leaving open the possibility of non singular superspinning solutions |7]. 

In a recent paper [8], the present authors showed that equilibrium configurations describing an extended compact object can 
in principle be obtained from gravitational collapse. The models we described correspond to a slowly evolving collapsing cloud 
which settles asymptotically to a static final configuration that is either regular or has a naked singularity at the center In this 
context, a naked singularity must be understood as a region of arbitrarily large density that is approached as comoving time goes 
to arbitrarily large values. We briefly examined the properties of accretion disks around the naked singularity solutions. There 
were significant physical differences compared to disks around black holes, and it followed that there could be astrophysical 
signatures that could distinguish black holes from naked singularities. 

In the above previous work, we considered the case of purely tangential pressure and vanishing radial pressure. Are the 
solutions we obtained a consequence of this extreme simplification, or are they representative of a generic class of solutions 
that survive even for more reasonable equations of state? We answer this question here at least partly by exploring gravitational 
collapse of perfect fluid objects with an isotropic pressure. Even in this more realistic case, we find that objects with either regular 
interiors or naked singularities form readily as a result of gravitational collapse. We explore the similarities and differences 
between these new solutions and the earlier solutions which had purely tangential pressure. 

Our main purpose in the present paper is to investigate whether the naked singularity models derived here and in our previous 
work can be observationally distinguished from black holes of the same mass. We therefore take a further step in this direction by 
calculating the spectral energy distributions of putative accretion disks. We show that important differences exist in the physical 
properties of accretion disks around naked singularities compared to those around black holes, which may help us distinguish 
black holes from naked singularities through observations of astrophysical objects. 

The structure of the paper is as follows: In section II we describe a procedure by which a static perfect fluid object with a 
Schwarzschild exterior metric can be obtained via gravitational collapse from regular initial conditions. In section III, we use 
the above procedure to obtain a toy model of a static final configuration with a naked singularity at the center We then describe 
the properties of accretion disks in this toy spacetime and compare these models to disks around a Schwarzschild black hole 
of the same mass. We also briefly discuss other density profiles of astrophysical interest that could be studied within the same 
framework. In the final section IV, we discuss possible applications to astrophysical observations. 

II. GRAVITATIONAL COLLAPSE 

Spherical collapse in general relativity can be described by a dynamical spacetime metric of the form 

7-)/2 

ds'^ ^ -e^^df + -^dr'^ + R^dVi^ , (1) 

where v, R and G are functions of the comoving coordinates t and r. For the perfect fluid case, the energy-momentum tensor is 
given by Tq = e, Tl — T2 — T^ — p. The Einstein equations then take the form 

TP 

(2) 



R^it ' 

F' 

(3) 



R'^R' ' 
P' 

e + p 



(4) 



G ^ 2^RG , (5) 



where (') denotes a derivative with respect to comoving radius r and (') denotes a derivative with respect to time t in the (r, t) 
representation. The function F, called the Misner-Sharp mass, describes the amount of matter enclosed within the sheU labeled 
by r at the time t, and is given by 

F = R{l-G + e-^''if). (6) 



Equations (|2|i-(|6]l give five relations among the six unknown functions p{r, t), e{r, t), v{r, t), G{r, t), F{r, t) and R{r, t). We 
thus have the freedom to specify one free function. An assumed equation of state relating pressure to energy density during 
the evolution of the system would fix this remaining freedom and give a closed set of equations. It could happen, however, 
that this approach leads to an analytically intractable problem. Moreover, we do not know if the collapsing matter will have 
the same unchanged equation of state as it evolves to higher and higher densities as the collapse progresses. There is in fact no 
astrophysical or mathematical requirement that the equation of state must be fixed as the collapse evolves in time. We therefore 
prefer here instead to choose the functional form of F{r, t), which corresponds implicitly to fixing the equation of state, which 
could however change with time and space coordinates. We choose a scenario such that we approximate some standard equation 
of state at early times, switching to some other, possibly exotic, kind of matter at later times. 

We use the scaling degree of freedom in the definition of R to fix the initial condition R{r, ti) — r, where ti is the initial 
time. To describe collapse we require R < 0, which guarantees that R decreases monotonically with respect to t. Hence, we 
may change coordinates from (r, t) to (r, R), thus in effect considering t = t{r, R). Correspondingly, we can view the original 
function F{r, t) as a function of r and R and write F = F{r, R). In the following, we use (, r) to denote a derivative with 
respect to r in the (r, R) representation, i.e., 

F' = F^r + FrR'. (7) 

The total mass of the system is not conserved during collapse, unless one requires the further condition that F{ri,, t) = const, 
where Vb corresponds to the boundary of the system. Therefore, we cannot in general match a collapse solution to an exterior 
Schwarzschild metric. However, matching to a generalized Vaidya spacetime at the boundary Rb{t) = R{ri,,t) is always 
possible [^. 

The procedure to solve the above system of Einstein equations is the following. We choose the free function F{r, R) globally 
and use equations (|2]l and (HJ to obtain p{r, R) and e(r, R). We then integrate equation (|4]i to obtain v, 

rr / 

z/(r,i?)=- / -^dP, (8) 
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and integrate equation (|5]l to obtain G, 

G(r,i?)=6(r)e2/"i5^'^-«. (9) 

The free function 6(r) results from the integral in equation (|9); it is related to the velocity profile in the collapsing cloud. The 
integral in equation (O again gives a free function of t, but this can be absorbed via a redefinition of the time coordinate. Once 
we have I'ir, R) and G{r, R), we can integrate the Misner-Sharp mass equation (|6}, written in the form 

^, (10) 
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to obtain tir, R), or equivalently R{r, t). 

Substituting equations (jSj and (|9]) in equation ([T]l, the metric of the collapsing spacetime takes the form 

ds^ = -e-^^o ^P^^dt^ + ^'l , ^ dr^ + R^dn" . (11) 

It may not always be possible to fully integrate the system of Einstein equations globally. However, this is not always needed, 
because by considering the behaviour of the functions involved, it is often possible to extract useful information about collapse 
and to integrate the solution at least in a neighborhood close to the center. 

Typically some restrictions are required in order for the collapse model to be considered physically viable. They are: 

1. Absence of shell crossing singularities, which arise from a breakdown of the coordinate system at locations where col- 
lapsing shells intersect. This requirement implies R' > 0. We note in this connection that there is always a neighborhood 
of the central line r = of the collapsing cloud which contains no shell-crossing singularities throughout the collapse 
evolution |10]. 

2. Energy conditions, of which the usual minimal requirement is the weak energy condition, viz., positivity of the energy 
density (e > 0) and of the sum of density and pressure (e + p > 0). This requires F'{r,R) > and F^r{r,R) > 0. 

3. Regularity at the center during collapse before the formation of the singularity. This includes the requirements that 
forces and pressure gradients vanish at the center and that the energy density has no cusps at r = 0. The corresponding 
requirements are F{r, R) ~ r^m{r, R) near r = 0, m'{0, R) — andp'(0, R) — 0. 



4. Absence of trapped surfaces at the initial time. This last requirement is given by positivity of [1 — {F/ R)]^^^ and 
ti-anslates to i?j^j^ < (e^^G)^^^ . 

Our aim now is to construct a dynamical collapse evolution such that the pressure eventually balances gravitational attraction 
and the collapsing object settles into a static configuration. Such a scenario is of course not always possible in gravitational 
collapse of a massive matter cloud, and there are matter configurations that can only collapse indefinitely without achieveing any 
possible equilibrium. Such is the case of collapse of a pressureless dust cloud flT], or similar systems where pressure is very 
insignificant. Another possibility in collapse is that the cloud bounces back after reaching a minimum radius. Typically, for each 
collapsing shell labeled by the comoving radius r, three different behaviours are possible: 

1. Complete collapse: 7? < at all times. The shell collapses to a central singularity. 

2. Bounce: 7? = at a certain time, with i? > at this time. The shell bounces back and re-expands. 

3. Equilibrium configuration: R = R = 0. Here the collapse slows down and the shell achieves a static configuration. 

Shells that achieve an equilibrium configuration in a certain sense mark the separation between the region that collapses indef- 
initely and the region that eventually bounces back fl2]. Both collapse and bounce would occur typically very rapidly, on a 
dynamical time scale which is proportional to the mass. Even in the case of supermassive compact objects, this time is much 
shorter than other typical time scales. 

It would appear that all three possibilities above represent generic behavior during gravitational collapse, depending on the 
masses and velocities of the collapsing shells and the physical scenario involved in collapse. For example, for a very massive 
star, an indefinite collapse would seem inevitable if the star cannot shed away enough of its mass in a very short collapse time 
scale to achieve any possible equilibrium. On the other hand, for a much larger mass scale such as galactic or yet larger scales, 
the collapse could proceed much more slowly and even arrive at an equilibrium to form a stable massive object. It is well known 
of course that gravitational collapse plays a key role in the formation of large scale structures in the universe. In many such 
cases, an evolving collapse would slow down, eventually to form a stable massive object. In such cases, scenarios such as the 
one considered here could be relevant. 

We note that in order to achieve an equilibrium object as the endstate of collapse or in order to have a useful quasi-equilibrium 
configuration, it is necessary that each collapsing shell in the cloud must individually achieve the condition R — R ~ 0. Such 
an object can be well approximated by a static configuration. The equation of motion ^ can be written in terms of R, for any 
fixed comoving radius r, in the form of an effective potential: 

V{r, R) = -if = -e"^" f^+G-l\ . (12) 

We already know that no static configuration is possible for the pressureless (dust) case, where V is negative at all times. When 
there is pressure, it is still possible for V to be negative at all times, giving continued collapse. However, other possibilities are 
also allowed. One generic possibility is that V has two distinct zeroes. In this case, the shell will bounce at a finite radius and 
re-expand. A third possibility is that the two zeroes coincide, corresponding to an extremum at y = 0. In this case, the shell will 
coast ever more slowly towards the radius corresponding to V^ = and will halt without bouncing. A global static configuration 
is then achieved only if each shell satisfies the condition that it has coincident zeros. It is not difficult to see that the velocity of 
the cloud during collapse is always non-zero and therefore such an equilibrium configuration can be achieved only in a limiting 
sense, as the comoving time t goes to infinity. 

The condition that the spacetime evolves towards an equilibrium configuration is thus 

R = R = 0, (13) 

which is equivalent to 

V ^Vr^O. (14) 

From equation (fT2l i we find 

The system achieves an equilibrium configuration if the solution of the equation of motion ( fTOl i. given by R{r, t), tends 
asymptotically to an equilibrium solution R^ (r) such that the conditions given in equations (fT3] l or (fl4l l are satisfied. Therefore, 
in order to have 

R{r,t) >Re{r), (16) 



we must choose the free function F{r, R) in the dynamical collapse scenario such that the quantities F, v, G tend to their 
respective equilibrium limits: 

F{r,R) -> F,(r)=F(r,i?e(r)), (17) 

u{r,R) ^ v,{r)^v{r,R,{r)), (18) 

G{r,R) -> G,{r)=G{r,R,{r)). (19) 

Imposing the conditions (fT4b . we thus obtain two equations which fix the behaviour of G and G^r at equihbrium: 

Ge(r) = 1 - ^ , (20) 



1 ^^ 




G,R{r,R,{T)) = 


Fe 

Ri 



(G,fl)e = G,r{t, R,{t)) = ^ ~ ^^^ , (21) 



where the velocity profile b{r) in equation (O has been absorbed into Ge{r). 
If we define the new radial coordinate p at equilibrium as 



we can rewrite the functions at equilibrium as 



p=Re{r), (22) 



Fe{r) = F(p), (23) 

Mr) - 4>{p), (24) 

Ge{r) - G(p)-1--. (25) 



P 
Then, from equations Q and (|4]i, two of the Einstein equations for a static source become 

e{p) = ^, (26) 

P^ 

P,P = -(e+p)0.p, (27) 

where (,p) denotes a derivative with respect to the new static radial coordinate p. The second equation is the well known 
Tolman-Oppenheimer-VoUcoff (TOV) equation. The third static Einstein equation, namely 

,.M.G(„-M ,28, 

P P-^ 

is obtained from equation dU by imposing the equilibrium condition and making use of equation (|5]l at equilibrium. The metric 
(fTTT i at equilibrium then becomes the familiar static spherically symmetric spacetime, 

ds^ = -e^^df + -^+ p^dn^, P<Pb^ Rein). (29) 

This interior metric is matched to a Schwarzschild vacuum exterior at the boundary pi, = Re{ri,). By matching gpp at the 
matching radius pi, and making use of equation ( 1251 ). we see that the total gravitational mass Mt of the interior is given by 

1-^ = GU) = 1-^M. (30) 

Pb pb 



Note that, in principle, the interior metric in equation (I29I I need not be regular at the center, as the eventual singularity is 
achieved as the result of collapse from regular initial data. A singularity at the center of the static final configuration is then 
interpreted as a region of arbitrarily high density that is achieved asymptotically as the comoving time t goes to arbitrarily large 
values. 

Typically, for a static perfect fluid source of the Schwarzschild geometry, various conditions could be required in order to 
ensure physical reasonability Ill3ll . Some of these are: the matter satisfies an energy condition, matching conditions with the 
exterior Schwarzschild geometry, vanishing of the pressure at the boundary, monotonic decrease of the energy density and 
pressure with increasing radius, sound speed within the cloud should be smaller than the speed of light. If the energy conditions 
are satisfied during collapse, they will be satisfied by the equilibrium configuration as well, i.e., the positivity of the energy 
density and sum of density and pressure at the origin follows from the same condition during collapse. Also, requiring only the 



weak energy condition allows the possibility of negative pressures either during the collapse phase or in the final equilibrium 
configuration. 

Since we have required here the static configuration to evolve from regular initial data, we may omit the condition that the 
final state must be necessarily regular at the center This leaves open the possibility that a central singularity might develop as 
the equilibrium is reached, where the singularity has to be understood in the sense explained above. Absence of trapped surfaces 
at the initial time ensures that, if a singularity develops, it will not be covered by a horizon lll4ll . In fact it can be shown that 
if trapped surfaces do form at a certain time as the cloud evolves, then the collapse cannot be halted and the whole cloud must 
collapse to a black hole or to a singularity where the first point of singularity is visible but the later portions of the singularity 
become covered in a black hole. We note that some of the above conditions, although desirable, may be neglected in special 
cases. Many interior solutions are available in the literature describing a static sphere of perfect fluid matching smoothly to a 
Schwarzschild exterior geometry (see [15] for a list of solutions or [16] for an algorithm to construct interior solutions). As 
discussed above, we can construct dynamically evolving collapse scenarios that lead asymptotically to the formation of such 
static configurations. If we are able to achieve a static configuration from collapse from regular initial data, and if we view a 
singularity as a region where density increases to arbitrarily high values, signaling a breakdown of the ability of general relativity 
to model the spacetime, we may then neglect the requirement of regularity at the center. 

III. A TOY MODEL OF A STATIC SPHERICALLY SYMMETRIC PERFECT FLUID INTERIOR 

In the following, we look for static interiors with a singularity at the origin. We wish to investigate the properties of such 
solutions and to establish whether such hypothetical naked singularity objects could in principle be distinguished in terms of 
observational signatures from black hole counter-parts of the same mass. 

We start with the most general static spherically symmetric metric written in the form given in equation (|29). As we have 
seen, the perfect fluid source Einstein equations give 






(31) 

-S£)l-I<£), ,32, 

P J P^ 

P,P = -(e+p)'/',p, (33) 

where we have absorbed the factor 47r into the definition of e and p, and defined F p = dF/dp. The third (Toknan-Oppenheimer- 
Volkoff) equation, combined with the second equation, gives 

^ [pp^+F{p)] 
'-^-^'^'^2p[p-Fip)Y ^''^ 

Since the above system of Einstein equations consists of three equations with four unknown, in order to close the system we can 
either specify an equation of state that relates p to e or use the freedom to choose arbitrarily one of the other functions IIITII . As 
discussed earlier, we opt for the latter approach and specify the form of the mass profile F(p), which describes the final mass 
distribution obtained from collapse. Once we specify F, equation ( |34| ) reduces to a first-order ordinary differential equation, 
which we can solve. 

We follow the same procedure that we used in the pure tangential pressure case ^ and obtain a solution of the form 

—^ = A/o = const, p< Pb- (35) 

P 

By equation dSOl l. this solution corresponds to a total mass Mt given by 

2Mt = F(pb) = Mopb. (36) 

Hence, to avoid a horizon, we require AIq < 1. Solving the first Einstein equation (l3Tb . the energy density becomes 

Mo 



P- 



,2 



(37) 



which clearly diverges as p — )• 0, indicating the presence of a strong curvature singularity at the center. 
The TOV equation (l34l i now becomes 



(e+p)2p _ (Mo+pp2)2 ^^^^ 



2(1 -Mo) 2p3(i_Mo)' 



which can be integrated by defining the auxiUary parameter 



Clearly, we require A G [0, 1), or equivalently AIq < 1/2 (for values of Mq > 1/2 similar considerations apply with A = 
■y/— (1 — 2Mo)/(l — Mq) but we do not consider this case here). In terms of A, the solution of equation dJSl l can be written as 

P ^ 2^^7 \ A^W' ) ' ^ ^ 

where A and B are arbitrary integration constants. Then, solving the remaining Einstein equation gives 

e^-^ = (V^^ - Bp^+^f. (41) 

This solution, among other similar interior static solutions, was first investigated by Tolman in 1939 for a specific choice of A 
01811 . The energy density of the solution may be rewritten as 

The existence of a strong curvature singularity at the center can be confirmed by an analysis of the Kretschmann scalar K near 
p = 0, which gives 



Q{p)^e"^^P\ (43) 



4p4Q4(2-A2)2 

which is clearly singular in the limit p — > 0. 

The static metric of the above solution takes the form 

ds^ = -{Ap^-^-Bp^+^fdt^ + {2-\^)dp^ +p'^dQ^, p<Pb. (44) 

It is matched at the boundary pi, to a vacuum Schwarzschild spacetime with total mass Mt- Since gpp in the interior is a constant, 
the matching will not in general be smooth, though it will be continuous. From the matching conditions for gpp and gu at the 
boundary we obtain, 

2(2 — A^) 

Pt = i: T^AfT, (45) 

(1- A^) 

A-l 

Bpl^ = A- ^'' (46) 

For all values of A £ [0, 1), the pressure is maximum at the center and decreases outwards, becoming zero at a finite radius pp 
given by 

p- ^ ^il^^ (47) 

^P B(l + A)2 ^ ^ 

Since we require the pressure to vanish at the boundary of the cloud, we obtain the further condition pb = pp, which, together 
with equation ( |46] |, fixes the two integration constants A and B: 

(1 + A)2p,^-i 

A = -^ 7=^, (48) 

4AV2^A2 

(1 - A)2p-^-i 

B = / ^'' (49) 

4A\/2^^A2 

The sound speed inside the cloud is given by c^ = dp/de. We find 

,^p AX^ABr^^ p ^ A'j^,-B^,p^^ 

^'^ e (l-A2)(^-Sr2A)2' e A-Bp^^ ' ^ ' 
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from which we see that Cs{0) = (1 — A)(l + A) < 1, as required, and Cs decreases as the radius p goes from to p^. We may 
also rewrite the above relations in the form of an equation of state, p = p{e), 



, , , l + A } \ 1-A2 } ^ ^ \ 1-A , 

PW = e^ ^7 k ^ -■ (51) 



m' 



In the limit A — > 0, corresponding to A/q -> 1/2, we obtain the most compact member of the above family of solutions. It 
has an equation of state p « e and hence sound speed c^ — > 1. Note, however, that at the outer boundary, p — > but e does 
not vanish, so the equation of state is not strictly isothermal. For this solution, the matching radius with the Schwarzschild 
exterior is at pt, ~ AMt, i.e., twice the Schwarzschild radius. In the case of pure tangential pressure, which we considered in 
our previous paper H, we found physically meaningful solutions down to pt — 3Mt, and even more compact solutions with 
exotic properties. This is one respect in which the perfect fluid model differs from the tangential pressure case. 

For A — 1/2, corresponding to Mq — 3/7, we obtain an equation of state 

P-ffl- ^ h (52) 

which approaches the radiation equation of state p = e/3 as p — > 0. For this model, pi, = (14/3)Mt, i.e., the object is a little 
less compact than the model with A — > 0. 

Finally, the case A -^ 1, Mo — > 0, corresponds to p^/Mt — > oo, and hence an infinitely large object. In this limit, our solution 
reduces to the classical Newtonian singular isothermal sphere solution. Correspondingly, A —> 1, B —^ 0, and the metric in 
equation ( |44] | reduces to the metric of flat space. 

Considering next the energy conditions, it is easy to see that A/q > implies positivity of the energy density and pressure. 
The pressure decreases from its maximum value at the center to zero at the boundary. From 

(1 -A)(1 + A) 1 f^ , p 
2 - A2 p 

we see that the weak energy condition is satisfied throughout the interior of the cloud. 



P= ' ^'\^ ^- 1 + -)' (53) 



A. Properties of circular orbits 

We wish to investigate basic observational properties of accretion disks orbiting in the above family of spacetimes. To this 
end, we study circular geodesies. Our assumption is that a test particle orbiting inside the "cloud" does not interact with the 
material of the cloud but merely feels its gravitational influence |32]. 

Our perfect fluid model has five parameters, Mt, A, pi,. A, B, and there are three matching conditions at the boundary p — pt,, 
viz., matching of gu and gpp with the exterior Schwarzschild metric and the condition p{pb) = 0. Thus, we are free to choose 
two of the five parameters. For convenience, in the following we choose the total gravitational mass of the cloud Mt to be equal 
to unity. This still leaves one free parameter, which we choose to be A. Once we pick a value of A, the boundary radius pb is 
given by equation (l45l l. and the coefficients A and B are given by equations ( l48T l and (l49T l. 

Given the spherical symmetry of the metric we can always choose the coordinate 6 such that the geodesic under consideration 
lies in the equatorial plane (9 — 7r/2). Time-like geodesies in this plane satisfy 

where 

H{p) = iAp'-'-Bp'+'f. (55) 

From the two killing vectors, ('', rj'^, associated with time-translational symmetry and rotational symmetry, we can calculate 
two conserved quantities, the energy per unit mass, E — gabC"'u'' — e^'^dt/dr, and the angular momentum per unit mass, 
L = gabV^u^ = p^dip/dr: 

dt 
E = H—, (56) 

dr 

L = p^^, (57) 

dr 



-1 — I — I — I — I — I — I — I — I — I — r 

2/ 




0.5 
log p 



FIG. 1: (a) Energy per unit mass E of circular orbits as a function of radius p for four perfect fluid models of unit mass (Mr = 1). The 
parameter A and the radius of the boundary Pi, of the four models are, respectively: (A, pb) = (1/100, 4.0002); (1/2, 14/3); (l/\/2, 6); 
(2/\/5, 12). The curves are labeled by their values of A. Each solution is matched to a Schwarzschild exterior al p — pt- (b) Close up of 
the matching region between the perfect fluid interior and the Schwarzschild exterior. The models with A = 1/100, 1/2, have pt lying inside 
the innermost stable circular orbit pisco = 6 of the Schwarzschild spacetime. For these two solutions, E decreases between pj, and pisco, 
indicating a zone of unstable circular orbits. 



For circular geodesies we must have dp/dr — 0. Therefore, we obtain 



p2 



2H^ 



2H-H^pp' 

H.pP 
2H - H^pp' 



(58) 
(59) 



Figure[T]shows the variation of the energy per unit mass E with radius p for a selection of models corresponding to A = 1/100, 
1/2, 1/V2, 2/a/5. The boundaries of the four models are at (see eq.l45t pb — 4.0002, 14/3, 6 and 12. In all the models, E goes 
to zero as p — > 0, with a power-law dependence: E ^ p^~^. At p = pb, each model is matched to an exterior Schwarzschild 
spacetime with mass Mt — 1. Panel (b) in Figure [T] shows the matching region. 

By construction, E is continuous across the matching boundary, but dE/dp is not. Note in particular that the vacuum 
Schwarzschild spacetime has an innermost stable circular orbit (ISCO) at pisco — 6. For p < pisco. E increases with 
decreasing p, which is one of the consequences of the absence of stable orbits. The model with A = 2/\/5 has its boundary at 
Pb — 12 and hence has stable circular orbits all the way from large radii down to p — >^ 0. So too does the model with A = l/\/2 
which has its matching radius at the ISCO, pb ~ Pisco- However, the other two models (A = 1/2, 1/100) have pb < pisco- 
Hence, these models have stable circular orbits inside pb and outside pisco = 6, but no stable orbits in between. 

Figure 12] shows analogous results for the angular momentum per unit mass L. Here, at small radii, L ^ p for all the models. 
As in the case of E, the angular momentum matches continuously across the boundary p — pb to the external Schwarzschild 
solution. As panel (b) shows, the two models with A = 1/100 and 1/2 have dL/dp < for a range of radii between pf, and 
Pisco- This corresponds to the region with unstable circular orbits. 

The angular velocity uj = d<j)/dt of particles on circular orbits is given by 



iU^ = 



9u,p 



H., 



This quantity scales as a; ^ p as p 



9<P<f>.,p 2p 
0. Figure|3]shows a plot of w vs p for the same four models as before. 



(60) 
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FIG. 2: (a) Angular momentum per unit mass L of circular orbits as a function of radius p for the same four perfect fluid models shown in 
Figure[T] The values of A are: 1/100, 1/2, l/\/2, 2/\/5 (the curves corresponding to the first and last solutions are labeled), (b) Close up of 
the matching region between the perfect fluid interior and the Schwarzschild exterior. Note, as in Figure [T] that the solutions with A = 1/100, 
1/2 have L decreasing between ph and pisco = 6. Circular orbits are unstable over this range of radius. 



B. Properties of accretion disks 



Gas in an accretion disk loses angular momentum as a result of viscosity and moves steadily inwards along a sequence of 
nearly circular orbits. Using just the properties of circular geodesies, and without needing to know the detailed properties of the 
viscous stress, it is possible to calculate the radiative flux emitted at each radius in the disk ll20ll . 

Since in our model both E and L tend to zero as p goes to zero, no energy or angular momentum is added to the central 
singularity by the gas in the accretion disk. The central singularity may be considered to be 'stable' in this sense. Indeed, since 
-E — > 0, all the mass energy of the accreting gas is converted to radiation and returned to infinity, i.e., the net radiative luminosity 
as measured at infinity satisfies Coo — ttk?, where m is the rate at which rest mass is accreted. Accretion disks around our 
model naked singularity solutions are thus perfect engines that convert mass into energy with 100% efficiency. 

From the behaviour of E and L of circular geodesies (Figs. [1] and |2}, we can distinguish two different regimes of accretion, 
depending on the value of A: 

• For A G [I/V2, 1), corresponding to pb £ [6Mt, +00), particles in the accretion disk follow circular geodesies of 
the Schwarzschild exterior until they reach the matching radius pi, at the outer edge of the cloud. Inside the cloud, the 
particles switch smoothly and continuously to the circular geodesies of the interior solution. Thus, the accretion disk 
extends without any break from arbitrarily large radii down to the singularity p ^ 0. 

• For A G (0, 1/a/2), corresponding to pf, G (4AiT, 6Mt), particles reach the ISCO of the exterior Schwarzschild spacetime 
at Pisco = 6. Inside this radius no stable circular orbits are allowed, so the gas in the disk plunges with constant 
E = Eisco and L = iisco until it crosses the boundary of the cloud at p^. Inside the cloud, circular geodesies are 
allowed again. The gas penetrates into the cloud until it reaches a radius pi at which the specific angular momentum of a 
local circular geodesic L is equal to Lisco- At this radius, the gas settles into a stable circular orbit and radiates away any 
excess energy. Further evolution than proceeds in the standard fashion, with the gas steadily moving to smaller radii until 
p —i' 0. The accretion disk is thus divided into two parts, one in the vacuum exterior over radii p > pisco and the other 
inside the cloud over radii p < pl- 

In the following, we focus on models belonging to the first regime, where we have a continuous disk extending with no 
break from large radius down to the singularity. Specifically, we consider models with A = l/V^, 2/\/5, which have matching 
boundaries at pf, = 6, 12, respectively. The radiative properties of accretion disks in these spacetimes may be calculated using 
the relations given in 1I20I1 . In the local frame of the accreting fluid, the radiative flux emitted by the disk (which is the energy 
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FIG. 3: (a) Angular velocity uj of circular orbits as a function of radius p for the same four perfect fluid models shown in Figure[T] The values 
of A are: 1/100, 1/2, l/\/2, 2/-\/5 (the curves corresponding to the first and last solutions are labeled), (b) Close up of the matching region 
between the perfect fluid interior and the Schwarzschild exterior. 



per unit area per unit time) is given by 

Hp) 



A-K^—g (E — UJ 



^ r {E-u:L)L.j,dp 



(61) 



where m is the rest mass accretion rate, assumed to be constant, pia is the radius of the inner edge of the accretion disk, which 
is zero for our singular cloud models, and g is the determinant of the metric of the three-sub-space {t, p, (/)), 



g{p) = -{2-X')p'H{p). 



(62) 



The solid curves in Figure |4la) show the variation of J-{p) vs p for the two chosen models. The flux diverges steeply as the 
gas approaches the center: J^ ^ p^^'^^'*''. This is not surprising, considering that the cloud is singular in this limit. Perhaps 
more surprising is the discontinuity in the flux at the boundary between the cloud and the external vacuum metric, as seen clearly 
in Figure |5b). While all the quantities E, L, to, g which are present in equation (l6ll are continuous across the boundary, the 
derivative w p is not. The discontinuity in w.p causes the jump in F as p crosses pt,. The dotted lines in the two panels correspond 
to an accretion disk around a Schwarzschild black hole. In this case, the inner edge of the disk is at pin = pisco = 6, and the 
flux cuts off at this radius. 

We note of course that the flux J^ is a local quantity measured in the frame of the fluid and is not directly observable. A more 
useful quantity is the luminosity Coo (energy per unit time) that reaches an observer at infinity. The differential of Coo with 
respect to the radius p can be computed from J^ by the following relation |I20']: 



dCr. 



d\n p 



A-np^/^EF. 



(63) 



Panels (c) and (d) in Figure|4]show dCoaldln p for the two models under consideration. We see that the luminosity behaves in 
a perfectly convergent fashion as p — > 0: dCoo/d\np ^ p^^^. By integrating this quantity over Inp, we can calculate the net 
luminosity Coo observed at infinity. We have confirmed that this is equal to m(? for the two singular models. That is the disk 
has 100% efficiency — it converts the entire rest mass energy of the accreting gas into radiation. The dotted lines in the panels 
indicate the very different behavior of a disk around a Schwarzschild black hole. Since such a disk is truncated at the ISCO, the 
luminosity is much less. In this case, we recover the standard result. Coo = (1 ^ -Bisco) "ic^ = 0.05719 mc^, i.e., an efficiency 
of5.719%. 

Is it possible to distinguish observationally whether a given accretion system is a black hole or one of the toy singular objects 
described in this paper? We have seen above that the accretion efficiencies are very different. However, the efficiency is not 
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FIG. 4: (a) The two solid lines show the variation of the radiative flux J" of an accretion disk, as measured by a local observer comoving with 
the fluid, for two perfect fluid singular models with A = l/\/2 (pb — 6) and A — 2/\/E (pb ~ 12). The dotted line corresponds to a disk 
around a Schwarzschild black hole, while the two dashed lines are for accretion disks in models with purely tangential pressure [8], which are 
discussed in section llil CI (b) Close-up of the flux profiles near the matching radius. Note the discontinuity in the flux, which is caused by the 
discontinuous behavior of uj^p in equation ( 161b . (c) Profile of the differential luminosity reaching an observer at infinity, dCoo/dlnp, for the 
same models, (d) Close-up of the region near the matching radius. 

easily determined via observations since there is no way to obtain an indepedent estimate of the mass accretion rate m. A more 
promising avenue is the spectral energy distribution of the disk radiation. 

Following standard practice, we assume that each local patch of the disk radiates as a blackbody. Defining a characteristic 
temperature T^, as follows (where we have included physical units), 



aTt = 



(64) 



47r(GAfT/c2)2' 
where a is the Stefan-Boltzmann constant, the local blackbody temperature of the radiation emitted at any radius p with flux T 
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FIG. 5: Spectral luminosity distribution of radiation from the accretion disk models in Figure |4] The dotted line corresponds to a disk around 
a Schwarzschild black hole, the solid lines correspond to disks around two perfect fluid naked singularity models discussed in this paper, and 
the dashed lines are for disks around two tangential pressure naked singularity models discussed in JSj]. 



(eq.|6T) is given by 



tMp) = [Hp)]'^'t, 



(65) 



This radiation is transformed by gravitational and Doppler redshifts by the time it reaches an observer at infinity and hence 
appears to have a different temperature. The transformation depends on the orientation of the observer with respect to the disk 
axis. To avoid getting into too much detail, we simply use a single characteristic redshift z, coiTesponding to an observer along 
the disk axis. 



l + z{p) = [-{gtt+uj^g^^)]-^/^, 
and assume that the radiation emitted at radius p has a temperature at infinity, independent of direction, given by 



(66) 



(67) 



It is then straightforward to convert the differential luminosity calculated in equation (|63]) into the spectral luminosity distribution 
'Ci/.oo observed at infinity. The result is 



r -i^ 



f dCo 



(1 + z)^{hv/kT:)^/F 
.^ \d\npj exp[(l + z)hi^/knJ^^/^] - 1 



dill p. 



(68) 



As a check, we have verified numerically that the integral over frequency of the spectral luminosity Lj^^oo obtained via the above 
relation is equal to mc^ (100% efficiency) for the singular models and equal to 0.05719mc^ for the Schwarzschild case (5.719% 
efficiency). 

Figure |5] shows spectra coiTesponding to the same models considered in Figure |4] For a given accretion rate m, the 
Schwarzschild black hole model gives a much lower luminosity than the naked singular models. However, as we discussed 
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earlier, this is not observationally testable. More interesting is the fact that the spectra have noticeably different shapes. At low 
frequencies, all the models have the same spectral shape vLi, ~ v^^^, which is the standard result for disk emission from large 
non-relativistic radii [21]. However, there are dramatic differences at high frequencies. 

The emission from a disk around a Schwarzschild black hole cuts off at the ISCO radius pisco = 6. Correspondingly, there 
is a certain maximum temperature for the emitted radiation, which causes the spectrum to cutoff abruptly at a frequency hv/kT^, 
somewhat below unity. The two perfect fluid singular models, on the other hand, behave very differently. In these models, 
the disk extends all the way down to p — > 0. Consequently, the temperature Tbb of the emitted radiation rises steadily and 
diverges as p — ^ 0. The emission from all the inner radii combines to produce a power-law spectrum at high frequencies [33]: 

As Figure |5] shows, the high-frequency power-law segment of the spectrum carries a substantial fraction (more than half) of 
the total emission from disks around our naked singularity models. The presence of this strong power-law spectrum is thus a 
characteristic feature of these models which can be used to distinguish them qualitatively from disks around black holes. Several 
well-known astrophysical black hole candidates are known to have spectra with a strong thermal cutoff, similar to that seen in 
the dotted line in Figure 15) ll22ll . These systems show very little power-law emission at high frequencies. In these cases at least 
we can state with some confidence that the central mass does not have a naked singularity of the sort discussed in this paper. 

C. Comparison of the perfect fluid and tangential pressure models 

In our previous paper H, we considered a relatively restricted model in which the fluid in the cloud has non-zero pressure 
only in the tangential direction. Because of vanishing radial pressure, matching pressure across the boundary radius pf, between 
the cloud and the external Schwarzschild spacetime was trivial, and hence the solution was more easily determined. For radii 
inside the cloud, the metric of the pure tangential pressure model takes the form, 

ds^ = -H{p)df+ \ +p^dn\ p<pi, (69) 

/ \ Afo/(l-Afo) 

H{p) = (1 - Mo) f -^ j , (70) 

Mo = ^, (71) 

Pb 

where Mt is the total gravitational mass of the cloud as measured in the vacuum exterior 

The close similarity of the tangential pressure model and the perfect fluid model described in the present paper is obvious, 
e.g., compare the above expressions with equations (l44l and (|36] |. The constant Mq in the tangential pressure model plays the 
role of A for the perfect fluid case, and in both cases this parameter determines the compactness of the cloud as measured by 
its dimensionless radius pi,/Mt- The limits Mq -> and A ^- 1 in the two models correspond to infinitely large and dilute 
clouds that are fully non-relativistic: pb/Mx -^ oo. In this limit, the two models are essentially identical. One minor difference 
between the models occurs in the opposite limit. For the perfect fluid cloud, the most compact configuration we find has A = 0, 
which corresponds to pb = 4Mt- In contrast, for the tangential pressure model, the most compact physically valid configuration 
has Mq = 2/3, which corresponds to pt, — 3Mt, i.e., a more compact object. 

Another difference between the two models is that the tangential pressure cloud has a particularly straightforward metric in 
which gtt — H{p) varies with radius as a simple power-law. Therefore, all quantities behave as power-laws and the analysis is 
easy. In the case of the perfect fluid model, the extra matching condition on the pressure at p = pi, results in the metric coefficient 
gtt involving two power-laws. The term involving the coefficient A in equation ( |44] | dominates as p ^- and behaves just like 
the lone power-law term in the tangential pressure cloud case. The second term involving B is required in order to satisfy the 
pressure boundary condition and plays a role only as p approaches pb. This term causes various quantities like E, L, uj, etc. to 
deviate from perfect power-law behavior as p — ?> pf, (see Figs. 1-4). 

As far as observables are concerned the two models behave quite similarly. In Figs. 4 and 5, the two perfect fluid models with 
A = I/V2 and 2/-\/5 have boundaries at pf, = 6 and 12, respectively, and their properties are shown by the solid lines. For 
comparison, the dashed lines show results for two tangential pressure models with Mq = 1/3 and 1/6 which have boundaries 
at the same radii, pi, — 6, 12. While the agreement between the two sets of models is not perfect, as is to be expected since 
the models are different, we see very good qualitative agreement. In particular, note that, just as in the perfect fluid models, the 
tangential pressure models too produce a strong power-law high energy tail in the spectra of their accretion disks (Fig. 5), which 
may be used to distinguish these models from disks around a Schwarzschild black hole (dotted line). 
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D. Comparison with the Newtonian singular isothermal sphere 

A simple and commonly used model in astrophysics is the singular isothermal sphere. This is a spherically symmetric self- 
gravitating object with an equation of state p = ec^, where p is the pressure, e is the density (this is usually called p but we use 
e here since we have already defined p to be the radius), and c^ is the isothermal sound speed. The model satisfies the condition 
of hydrostatic equilibrium under the action of Newtonian gravity. The solution is 

^^^=2^^ P^^^-2^^ M(P) = ^, (72) 

where M{p) is the mass interior to radius p. The variations of e, p and M with p are very reminiscent of the variation of e, p and 
F in the perfect fluid relativistic cloud model described earlier. 

One deficiency of the basic singular isothermal sphere model is that it extends to infinite radius, where the mass is infinite. In 
order to obtain an object with a finite radius, one needs to change the equation of state such that the pressure goes to zero at a 
finite density. This is easily arranged as follows: 

c o c 

e(p) = ^ttG 2 "" ^'^'" ^^ = 2ttG ^ ' ^^^^ "" ^"'-'^ ~ '^''^^ ^'^^^^ " 47ref,rgp, p < pt, (73) 

where pb is the radius of the boundary and £{, is the density at that radius. 

Circular orbits inside a singular isothermal sphere behave very simply, with velocity and angular momentum given by 

fcirc = \/crp = const, L = V2csP, P < Pb- (74) 

The constancy of Vdic is the chief attraction of the singular isothermal sphere. It provides a simple way of reproducing the 
observed flat rotation curves of galaxies. The scaling of L with p in the singular isothermal sphere is exactly the same as the 
variation of L with p in the perfect fluid relativistic cloud. At first sight it appears that the flat rotation curve property is not 
reproduced in the relativistic model. For example, ujp scales as p^^'*'. However, note that uj is defined as d4)/dt, where t is the 
time measured at infinity. The appropriate quantity to consider is d(j)/dt\oc, where tioc is time measured by a local ZAMO (Zero 
Angular Momentum Observer) at radius p. The two times are related by a factor of [H{p)] ^^^ which varies as Ap^^^ in the limit 
p — > 0. Thus, for ZAMOs, the rotation velocity is indeed independent of radius in the deep interior of the cloud. (Close to the 
boundary, there are small deviations because of the term involving B in H{p)). 

Consider next the energy of circular orbits in the singular isothermal sphere model. If we include the rest mass energy and 
add to it the orbital kinetic energy and the potential energy, then a Newtonian calculation gives for the energy of a particle of 
unit mass 

^circ = 1 - ^ ("l + 2 In ^ V p<pt,. (75) 

The weak logarithmic divergence at small radii is rarely a problem since we generally have Cg ^ c. Nevertheless, the presence 
of the logarithm implies that, in principle, at a sufficiently small radius, the Newtonian model predicts a negative total energy for 
the particle. This is clearly unphysical. 

It is reassuring that the logarithmic divergence is not present in our relativistic cloud models. Both the pefect fluid model 
and the tangential pressure model have E varying as a power-law with radius: E ^ p^^^ for the perfect fluid case and E ^ 
pMo/2{i-Mo) fQj. jjjg tangential pressure case. In both cases, E asymptotes precisely to zero as p — > and does not go negative 
anywhere. Effectively, the relativistic models, being more self-consistent, regularize the logarithm of the Newtonian model by 
replacing it with a power-law. The index of the power-law is nearly in the Newtonian limit but becomes as large as unity for the 
maximally compact configuration, viz., A = 1 and Mq — 2/3 for the perfect fluid and tangential pressure models, respectively. 

E. Other models 

As we have seen above, the perfect fluid relativistic model described in this paper is equivalent to the Newtonian singular 
isothermal sphere model used in astrophysics. The relativistic generalization was obtained by assuming that the energy density 
e(/5) in the relativistic model has the same functional form as e(p) in the Newtonian model (compare eqs.[37landl73Tl. and then 
solving for the pressure p and the spacetime metric. This procedure can be followed with any other trial model of e(p) of interest. 

One simple example is to consider a superposition of the toy perfect fluid model described in this paper with the well known 
constant density Schwarzschild interior. The density profile then takes the form 

e=^ + Mi, (76) 
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which corresponds to 

F{p) = Mop + AfipVS. (77) 

The density in this model approaches the singular interior of our perfect fluid model as p — > but resembles a constant density 
interior with e ~ Mi as p ^- pb. It turns out that the TOV equation can be explicitly integrated in this case, though the expression 
forp is fairly complicated and involves hypergeometric functions. 

Other examples of more interest to astrophysics could be similarly considered. One natural generalization of the singular 
isothermal sphere is the Jaffe density profile 112 311 . 

- ^'°'^° (78) 



where Mq is a constant and po describes the characteristic radius of the object. This density profile corresponds to 

F{P) = j^. (79) 

(Po + P) 

The Jaffe model behaves just like the singular isothermal sphere as p — > 0, yet it has a finite total mass given by Mq and hence 
does not need to be artifically truncated as we had to in the case of the singular isothermal sphere. The Jaffe model is a special 
case of a more general class of models, the Dehnen density profile ll24ll . 

^ (3-7)MoPo (30) 



p^(p + P0)4-7' 

The Jaffe profile corresponds to 7 = 2, while the case with 7 = 1 is known as the Hernquist model |25]. The Dehnen profile 
implies a mass function 

F{p) = Afo — ^ . (81) 

\Pa+ Pj 



(82) 



Finally, we could also consider the Navarro-Frenk- White (NFW) profile Il26l] . which is given by 

p(po + pY ' 



where again A/q is a constant. This model has a logarithmically diverging mass as p — > 00, and is thus a little less attractive. 

It is not easy to solve the TOV equations for the pressure p and the spacetime metric of any of the above models analytically. 
However, numerical solutions are easily obtained. Note that all these popular models technically have naked singularities at their 
centers. Obtaining relativistic generalizations along the lines followed in this paper would be worthwhile. 

IV. CONCLUDING REMARKS 

In the present paper we considered a self-bound spherical cloud of perfect fluid and derived static non-vacuum solutions of the 
Einstein equations which posseses a central naked singularity. We showed that these solutions could be obtained asymptotically 
as the final result of the slow collapse of a massive matter cloud. The solutions described here closely parallel those we obtained 
in ||8|] for a fluid with pure tangential pressure. 

We studied the properties of steady thermal accretion disks in our naked singularity spacetimes. Focusing on those models 
that have a disk extending continuously from a large radius p down to p — > 0, i.e., models with A > 1/a/2 for the perfect fluid 
case and Mq < 1/3 for the tangential pressure case, we showed that accretion disk spectra would consist of a multi-temperature 
blackbody at low frequencies joining smoothly to a power-law at high frequencies. Notably, the disk luminosity would be 
dominated by the high-energy tail, which is a characteristic feature of these models. The spectrum of an equivalent accretion 
disk around a black hole would have only the low-frequency multi-temperature blackbody component and would be missing the 
high frequency tail (or at best have a weak tail). 

Accretion disks around astrophysical black hole candidates do show power-law tails, but these are usually interpreted as 
coronal emission from hot gas above the disk. In those cases where the disk emission is definitely thermal, e.g., in the so-called 
Thermal State of accreting stellar-mass black holes [27], the spectrum is invariably dominated by the low-frequency multi- 
temperature emission and the power-law tail tends to be quite weak. A number of stellar-mass black holes have been observed 
in the Thermal State |28]. In the case of these systems at least we may conclude that the central compact objects are not naked 
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singularities of the type discussed in the present paper; the objects are presumably true black holes. However, it is not possible 
to state anything with certainty in the case of other black hole candidates that have not been observed in the Thermal State. 

The existence of horizons in stellar mass astrophysical sources has been a matter of discussion in recent years (see for example 
II29II ). Some observations suggest that the departure of such objects from black holes must be very small |30]. All the same, 
when it comes to more massive sources, excess luminosity from Ultra Luminous X-ray Sources cannot at present be explained 
fully by the usual black hole accretion disk models and seem to require the existence of as yet undiscovered intermediate mass 
black holes iH. 

Since there is a strong thermal cutoff at high frequencies present in the spectrum of the black hole models while it is absent in 
the perfect fluid naked singularity models studied here, this strenghtens the argument that the astrophysical sources that exhibit 
the same behaviour would be black holes. On the other hand, ultra luminous soruces exist in the universe and at present it is 
not clear if it will be possible to fit all the observations within the black hole paradigm. In such a scenario, the comparison and 
models such as discussed here may be useful to understand better and analyze. 

For the toy models presented here we found that the radiant energy flux and the spectral energy distribution are much greater 
as compared to a black hole of the same mass, and therefore we obtained indications that if some sources of similar kind do exist 
in the universe it might be possible to distinguish them observationally from black holes. 

We also found here that the qualitative features described in |8] are preserved when we consider perfect fluid sources, instead 
of sources sustained only by tangential stresses. We can therefore conjecture that the increased flux and luminosity may be a 
generic feature of any source with a singularity at the center where the accretion disk can in principle extend until r = 0. As it 
is known, perfect fluids are important in the context of astrophysics as they can be used to model many sources and objects of 
astrophysical relevance. Of course, different matter models will imply different luminosity spectrums that might or might not 
eventually be distinguished from one another 

Finally, we note that these models require a certain fine tuning since all the collapsing shells must have the right velocity 
so that the effective potential leads asymptotically to R ^ R = 0. Nevertheless the class of static final configurations can be 
understood in the sense of an idealized approximation for a slowly evolving cloud with small velocities. In this sense these 
scenarios appear to be generic as they approximate a wide variety of slowly evolving matter clouds. 
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